function f = elliptic_f(x,y,caseNum,varargin)
global alpha theta;
        f=(x.^2+y.^2).^(3/alpha).*sin((2*atan2(y+eps,x+eps)+theta)/alpha).^3;
   